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EXECUTIVE  SUMMARY 

The  primary  goal  of  this  research  was  to  develop  physically  based  damage  models, 
computational  approaches  based  on  shear  deformation  plate  theories  and  layerwise  theory,  a 
finite  element  formulation  for  the  sinalysis  of  thin  and  thick  laminated  composite  structures 
with  embedded  sensors/actuators.  During  the  period  of  this  project  the  following  tasks  have 
been  completed:  (1)  literature  search  to  assess  damage  and  failures  models;  (2)  theoretical 
formulation  of  the  layerwise  models  for  smart  composite  structures,  (3)  formulation  of  damage 
models,  (4)  finite  element  analysis  of  laminated  adaptive  plate  structures  using  the  layerwise 
theory,  and  (5)  finite  element  analysis  of  laminated  adaptive  plate  structures  using  the  sin^e- 
layer  theories  (classical,  first-order,  and  third-order  plate  theories).  Geometric  nonlinearities 
are  accounted  for  and  parametric  effects  of  lamination  schemes,  material  properties,  and 
boundary  conditions  were  investigated. 


TECHNICAL  DISCUSSION 


1.  Damage  Models 
1.1  Introduction 

Fiber-reinforced  laminated  composites  typically  exhibit  several  different  forms  of  progressive 
damage  including  fiber  breakage,  fiber  buckling,  fiber/matrix  debonding,  matrix  cracking, 
individual  ply  rupture  and  separation  between  adjacent  plies  (i.e.  delamination.)  These  modes 
of  damage  are  often  initiated  at  loads  that  are  far  below  the  ultimate  failure  load  of  a  composite 
material.  Any  damage  that  is  initiated  tends  to  progress  as  the  loading  is  further  increased  or 
cycled,  and  thus  effectively  degrades  the  stiffness  of  the  material,  causing  significant  load 
redistribution  before  ultimate  failure  occurs.  To  capture  this  complex  behavior  in  a  finite 
element  simulation,  two  different  approaches  have  emerged  for  the  development  of  damage- 
dependent  constitutive  relations:  a  phenomenological  approach  (known  as  continuum  damage 
mechanics,  or  CDM),  and  a  micromechanics  approach.  In  practice,  the  difference  between 
these  two  approaches  is  often  blurred;  however,  in  the  extreme  case,  we  can  differentiate 
between  the  two  approaches  as  follows. 

•  The  CDM  approach  attempts  to  predict  damage  accumulation  and  stiffness  reduction  in 
the  equivalent  homogenized  version  of  the  heterogeneous  material. 

•  The  micromechanics  approach  attempts  to  predict  damage  accumulation  and  stiffness 
reduction  in  the  heterogeneous  material  before  homogenization. 
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Therefore  we  can  say  that  the  difference  between  the  two  approaches  is  simply  a  matter  of  the 
timing  of  the  homogenization  process.  The  micromechanics  approach  has  the  potential  for 
greater  accuracy;  however,  it  typically  involves  significantly  more  computational  effort  and  also 
requires  more  experimentally  derived  material  data  (much  of  which  is  often  unavailable).  Given 
these  obstacles,  "from  a  practical  point  of  view,  it  is  difficult  to  imagine  that  a  purely 
micromechanical  theory,  with  all  of  its  complexities,  will  ever  replace  a  properly  formulated 
phenomenological  theory  as  a  design  tool,"  (Krajcinovic  [1,2]).  Perhaps  in  its  most  useful 
sense,  the  micromechanics-based  damage  model  can  be  utilized  as  a  basis  on  which  a 
simplified  COM  model  can  be  built  with  a  reasonable  chance  of  replicating  experimentally 
observed  trends  over  a  wide  spectrum  of  materials  and  responses  (Krajcinovic  [1,2]).  This  is 
essentially  the  same  philosophy  adopted  in  the  current  proposal  for  a  general-purpose  finite 
element  framework  capable  of  accurately  and  efficiently  simulating  the  progressive  damage  of 
smart  composites. 

1.2  Literature  Review 

A  detailed  literature  review  of  the  existing  damage  and  failure  models  was  included  in  the  first 
interim  progress  report.  Here  we  briefly  discuss  the  main  developments  of  that  review. 

Chang  and  Chang  [3]  gave  a  progressive  damage  model  for  laminated  composites  containing 
stress  concentrations  and  subjected  to  tensile  loading.  The  aim  was  to  assess  damage  and 
predict  ultimate  tensile  strength.  Stress  and  strain  analysis  was  based  on  classical  laminate 
theory,  with  material  non-linearity  also  being  considered.  For  matrix  cracking  failure,  they 
proposed  a  matrix  failure  criterion.  Fiber-matrix  shearing  and  fiber  breakage  were  predicted  by 
a  modified  Yamada-Sun  failure  criterion.  A  property  degradation  model  was  also  proposed. 
Results  obtained  from  a  nonlinear  finite  element  model  for  laminates  containing  a  circular  hole 
were  compared  with  experimental  data. 

Aboudi  [4]  presented  a  model  in  which  the  effect  of  damage  due  to  imperfect  bonding  between 
constituents  of  composite  materials  was  incorporated.  The  interface  de-cohesion  was  described 
by  two  parameters  which  completely  determine  the  degree  of  adhesion  at  the  interfaces  in  the 
normal  and  tangential  directions.  Perfect  bonding,  perfect  lubrication  and  complete  debonding 
were  obtained  as  special  cases.  The  model  also  predicts  the  overall  modulii  and  overall 
coefficients  of  thermal  expansion  of  composites  in  the  presence  of  imperfect  bonding.  Ladeveze 
and  Le  Dantec  [5]  presented  a  model  of  the  mechanical  behavior  of  am  elementairy  ply  in  a 
fibrous  composite  laminate.  Continuum  damage  mechanics  theory  was  used  to  describe  matrix 
microcracking  amd  fiber-matrix  debonding.  To  model  the  anelastic  strains  induced  by  damage, 
a  plasticity  model  was  proposed.  The  model  takes  into  account  the  differences  between  tension 
amd  compression  in  the  fiber  direction.  The  model  leads  to  a  laminate  failure  criterion. 

Homogenization  approaches  to  damage  have  also  been  developed.  These  often  involved  a  two- 
step  process  in  two  different  ways:  a)  homogenization  of  the  stationary  (non-evolving) 
microstructure,  i.e.  fibers  and  matrix,  followed  by  incorporation  of  the  evolving  damage  as 
internal  variables,  and  b)  incorporation  of  damage  in  each  phase  followed  by  homogenization  of 
the  phases.  Fish  et  al.  [6]  choose  to  carry  out  the  two  steps  simultaneously  by  using  an 
extended  homogenization  procedure,  wherein  the  field  variables  and  a  damage  variable  are 
given  asymptotic  expansions.  The  micro  level  total  displacement  is  assumed  to  have 
contributions  from  the  elastic  strain  and  a  damage  induced  strain,  d,  which  is  related  to  the 
damage  variable  (O.  The  microstructure  of  a  composite  material  is  assumed  to  be  locally 
periodic  and  the  influence  of  this  periodicity  scale  is  assumed  to  be  the  same  for  both  strains. 
To  avoid  a  computational  stability  problem  observed  previously  in  local  approaches.  Fish  et  al. 
defined  a  non-local  damage  variable  that  is  a  weighted  average  over  the  two  phases  in  a 
characteristic  volume.  Finally,  they  expressed  damage  evolution  in  each  phase  separately  in 
terms  of  the  damage  energy  release  rates.  Another  example  of  the  homogenization  approach  is 
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that  of  Ghosh  and  his  colleagues  [7-9].  These  authors  do  not  use  the  concept  of  damage  as  an 
internal  variable;  instead,  they  determine  the  damage  events  explicitly  by  local  stress  analysis 
and  failure  criteria.  They  also  allow  for  non-uniformity  in  the  distribution  of  one  phase 
(particles  or  fibers)  into  another  (matrix),  and  use  Voronoi  cells  to  represent  this  heterogeneity. 
The  micro  scale  model  is  analyzed  using  the  effective  stradn  from  the  global  (structure)  scale 
results  and  the  effective  (average)  modulus  tensor  is  evaluated.  When  damage  evolves,  the  two- 
scale  analysis  must  be  conducted  iteratively. 

One  of  primary  difficulties  that  one  faces  in  attempting  to  model  progressive  damage  in  fiber- 
reinforced  laminated  composites  is  that  many  of  the  damage  modes  (e.g.  fiber  breakage,  fiber 
buckling,  fiber/matrix  debonding,  matrix  cracking)  occur  on  geometric  scales  that  are 
impractical  to  reach  with  direct  mesh  refinement.  However,  these  forms  of  damage  have  a  very 
significant  effect  on  the  stiffness  of  the  heterogeneous  material  that  constitutes  the  typical 
finite  element.  On  the  other  hand,  some  forms  of  damage  (e.g.  ply  rupture  and  delamination) 
occur  on  geometric  scales  that  might  be  smaller  or  larger  than  the  t3TJical  finite  element, 
depending  on  the  level  of  mesh  discretization.  Furthermore,  the  criteria  used  to  evaluate  the 
onset  and  progression  of  damage  invariably  involves  the  stress,  strain  and/or  energy  of  the 
materieil  and  these  results  are  known  to  exhibit  a  strong  dependence  on  mesh  descritization. 
Therefore,  in  attempting  to  develop  a  computational  model  that  incorporates  all  these  various 
forms  of  damage,  we  must  distinguish  between  macroscopic  damage  and  microscopic  damage. 
Within  the  context  of  modeling  and  simulation,  we  will  distinguish  between  these  two 
categories  of  damage  on  the  basis  of  whether  the  damage  entity  is  smaller  or  larger  than  the 
computational  reference  scale  (i.e.  the  geometric  scale  of  the  individual  finite  element.)  Thus, 
our  choice  of  using  the  geometric  scale  of  the  individual  finite  element  as  the  reference  for 
distinguishing  microscopic  damage  from  macroscopic  damage  leads  to  the  following 
interpretations: 

•  Modeling  of  microscopic  damage  does  not  explicitly  require  modification  of  the  basic 
assumed  kinematics  of  the  model,  but  does  require  modification  of  the  materiEil  properties 
to  account  for  the  damaged  or  weakened  state  of  the  material.  Typically,  microscopic 
damage  is  handled  by  using  a  homogenized  description  of  the  distribution,  density  and 
orientation  of  the  microcracks  and  microvoids  that  occur  in  the  immediate  neighborhood  of 
a  Gaussian  integration  point  within  a  finite  element.  For  finite  elements  that  are 
significantly  larger  than  the  fiber/matrix  RVE,  but  not  necessarily  lairger  than  the  thickness 
of  a  single  lamina,  microscopic  damage  typically  includes  fiber  breakage,  fiber  buckling, 
fiber/matrix  debonding  and  matrix  cracking. 

•  Modeling  of  macroscopic  damage  explicitly  requires  modification  of  the  basic  assumed 
kinematics  of  the  model.  This  category  of  damage  includes  cracks  or  voids  that  are  of  the 
same  order  or  lairger  than  the  representative  finite  element  size.  For  finite  elements  that  are 
smaller  than  the  thickness  of  a  single  laimina,  macroscopic  damage  might  include  ply 
rupture  and  delamination.  However,  for  very  large  finite  elements  whose  smallest 
dimension  is  equivalent  to  the  laminate  thickness,  even  these  forms  of  damage  might  be 
treated  as  microscopic  damage.  Of  course,  the  appearance  of  macroscopic  damage  is 
necessarily  preceded  by  cumulative  microscopic  damage  that  coalesces  into  macroscopic 
damage.  Therefore,  a  simulation  that  is  capable  of  predicting  the  appearance  and 
subsequent  propagation  of  macroscopic  damage  must  necessarily  incorporate  both  a 
mechanism  to  alter  the  material  properties  based  oh  homogenized  microscopic  damage  and 
a  mechanism  to  alter  the  assumed  kinematics  based  on  macroscopic  cracks  and  voids. 

From  the  above  description,  it  is  obvious  that  the  classification  of  a  particular  form  of  Hamagp 
as  microscopic  damage  or  macroscopic  damage  depends  upon  the  specific  type  and  quality  of 
computational  results  that  we  seek.  This  is  the  most  influential  factor  in  our  choice  of  a 
particular  mathematical  model  and  a  particular  finite  element  size.  Thus  it  is  imperative  that  a 
general-purpose  computational  tool  should  be  flexible  in  its  treatment  of  damage  modeling  at 
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different  geometric  scales.  Ideally,  the  classification  of  different  physical  modes  of  damage  as 
microscopic  or  macroscopic  should  be  automatically  based  on  the  characteristic  element  length 
scale  that  is  appropriate  for  the  particular  results  sought.  Thus  for  one  type  of  problem,  the 
simulation  code  might  handle  delamination  as  a  macroscopic  form  of  damage,  while  for 
another  type  of  problem  the  code  might  handle  delamination  as  a  microscopic  form  of  damage. 
While  the  complete  satisfaction  of  this  goal  is  impractical  given  the  current  state  of  simulation 
technology,  the  software  development  proposed  in  Phase  II  of  this  effort  represents  a  significant 
step  in  this  direction. 


1.3  Microscopic  and  macroscopic  damage  modeling 

The  computational  model  for  simulating  progressive  damage  in  smart  composites  incorporates 
the  following  mechanisms  for  modeling  the  initiation  and  propagation  of  damage  induced  by 
service  loads. 

1.  Microscopic  damage  is  modeled  within  the  context  of  Continuum  Damage  Mechanics,  or 
CDM,  where  the  material  properties  are  modified  to  account  for  the  weakened  state  of  the 
material  due  to  the  presence  of  distributed  microcracks  and  microvoids  that  occur  as  a 
result  of  stress  concentrations  within  the  heterogeneous  laminate.  With  the  exception  of 
pre-existing  delaminations,  all  forms  of  damage  are  initially  modeled  as  microscopic 
damage  using  continuum  damage  mechanics,  regardless  of  the  size  of  a  particular  finite 
element.  This  method  is  appropriate  up  to  the  point  where  the  microscopic  damage 
becomes  severe  enough  to  cause  coalescence  into  macroscopic  damage  (i.e.  delamination  or 
ply  rupture)  whose  dimension  is  approximately  equal  to  the  dimension  of  the  finite  element. 
The  successful  implementation  of  the  modeling  of  microscopic  damage  via  continuum 
damage  mechanics  requires  the  following  components: 

A.  Damage  Variable  -  Choose  an  appropriate  tensorial  form  for  the  internal  variable  that 
will  be  used  to  characterize  the  microscopic  damage.  This  internal  variable,  known  as 
the  damage  tensor,  is  responsible  for  describing  the  distribution,  density,  and 
orientation  of  microcracks  and  microvoids.  Various  forms  of  the  damage  tensor  have 
been  proposed  in  the  literature,  for  example,  scalars,  vectors,  2"<*  order  tensors,  and  4* 
order  tensors  (Skrzypek  and  Ganczarski  [10]).  Depending  on  the  anticipated 
complexity  of  the  damage  process,  we  should  choose  the  simplest  tensor  that  is  capable 
of  accurately  describing  the  distribution,  density,  and  orientation  of  microcracks  and 
microvoids  for  the  intended  type  of  problem.  For  fiber-reinforced  composite  laminae, 
the  simplest  form  of  the  damage  tensor  that  is  capable  of  accurately  describing 
microscopic  damage  is  a  ssnnmetric  2"“*  order  tensor  whose  principal  directions  are 
assumed  to  coincide  with  the  principal  material  directions,  i.e.  orthotropic  damage 
(Barbero  and  De  Vivo  [11]). 

B.  Stiffness  Reduction  Scheme  -  Define  the  effective  homogenized  material  properties  in 
terms  of  the  current  level  of  the  damage  variable,  i.e.  we  must  have  a  consistent  means 
to  modify  the  original  homogenized  material  properties  to  reflect  the  wesikened  state  of 
the  material  due  to  the  presence  of  distributed  microcracks  and  microvoids.  This  is 
accomplished  by  using  various  equivalence  principles,  e.g.  the  principle  of  elastic  energy 
equivalence. 

C.  Damage  Surface  -  Define  a  damage  threshold  criterion,  i.e.  a  multidimensional  surface 
in  the  space  of  state  variables  that  must  be  exceeded  before  any  damage  can  be 
initiated,  or  before  existing  damage  can  progress.  Thus  the  damage  surface  serves  to 
separate  non-damaging  behavior  from  damage-inducing  behavior. 


5 


D.  Damage  Evolution  -  Define  the  damage  evolution  equation.  Assuming  that  the  damage 
surface  has  been  reached,  we  must  express  the  rate  of  change  of  the  damage  tensor  as 
a  function  of  the  remaining  state  variables. 

2.  Macroscopic  damage  in  the  form  of  delamination  is  modeled  by  hierarchically  enriching  the 
kinematics  of  the  model  to  permit  general  relative  motion  between  the  newly  created 
interlaminar  surfaces.  Thus  the  appearance  of  macroscopic  damage  involves  an  increase  in 
the  number  of  degrees  of  freedom  of  the  model.  The  successful  implementation  of  the 
modeling  of  macroscopic  damage  requires  the  following  components: 

A.  Kinematic  Enhancement  -  Incorporate  a  hierarchical  description  of  the  displacement 
field  that  permits  degrees  of  freedom  to  be  added  to  the  existing  model  to  account  for 
the  relative  motion  of  newly  created  interlaminar  surfaces  due  to  delamination.  The  use 
of  a  hierarchical  displacement  field  facilitates  automation  of  the  process  since  the 
existing  mesh  does  require  alteration;  we  simply  add  incremental  displacement  degrees 
of  freedom  to  the  existing  model  (Robbins  and  Reddy  [5,6]). 

B.  Delamination  Criterion  -  Define  a  criterion,  based  on  the  homogenized,  microscopic 
damage  variable  to  indicate  initiation  and  further  propagation  of  delamination. 

A  detailed  description  of  modeling  of  microscopic  and  macroscopic  damage  will  be  given  in  the 
following  sections. 

1.4  Continuum  damage  mechanics  at  the  meso-scale  (Le.  ply  Level) 

In  the  sequel,  the  use  of  the  term  damage  is  understood  to  mean  microscopic  damage  unless 
explicitly  stated  otherwise.  Thus  the  term  damage  is  intended  to  indicate  the  presence  of 
microcracks  and  microvoids  that  are  small  compared  to  the  geometric  dimensions  of  a  typical 
finite  element  and  can  thus  be  handled  to  sufficient  accuracy  via  continuum  damage 
mechanics,  or  CDM.  In  most  practical  analyses  of  composite  materials,  the  typical  dimensions 
of  a  finite  element  do  not  approach  the  geometric  scale  of  the  fiber/matrix  RVE.  Thus  it  can  be 
safely  assumed  that  cracks  and  voids  that  are  small  compared  to  the  thickness  of  a  single 
fiber-reinforced  composite  lamina  will  be  treated  as  microscopic  damage  and  handled  via  CDM. 
Thus  damage  in  the  form  of  fiber  breakage,  fiber/matrix  debonding,  and  matrix  cracking  will 
be  treated  via  homogenization  at  the  meso-scale  (or  ply  level).  This  is  analogous  to  the  process 
used  to  define  equivalent  anisotropic  material  properties  at  the  meso-scale  (or  ply  level)  by 
homogenization  of  the  heterogeneous  material  constituents.  In  contrast,  the  modeling  of 
delamination  most  often  requires  treatment  as  a  form  of  macroscopic  damage  that  cannot  be 
homogenized. 

The  meso-scale  continuum  damage  mechanics  model  developed  in  this  study  is  capable  of 
utilizing  the  full  3-D  stress  and  strain  fields  evaluated  at  the  ply  level.  This  model  represents 
an  extension  of  the  2-D  continuum  damage  mechanics  model  of  Barbero  and  DeVivo  [11], 
where  each  individual  lamina  was  assumed  to  be  in  a  state  of  2-D  plane  stress.  In  the  present 
CDM  approach,  the  fiber-reinforced  composite  material  that  constitutes  an  individual  lamina  is 
replaced  with  a  non-heterogeneous  orthotropic  material.  Prior  to  the  appearance  of  any 
damage,  the  equivalent  orthotropic  material  may  be  considered  homogeneous  within  the 
lamina.  However,  once  damage  initiates,  the  orthotropic  material  properties  vary  with  respect 
to  position  within  the  lamina  since  damage  occurs  locally  within  the  lamina  depending  on  the 
local  state  of  the  material. 

1.5  Description  of  damage  (kinematics) 
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A  symmetric  2“d  order  damage  tensor  D  =  Dyetej  is  chosen  to  describe  the  distribution,  density 
and  orientation  of  microcracks  and  microdefects  within  the  homogenized  lamina  When  the 
damage  tensor  D  is  expressed  in  an  arbitrary  coordinate  system,  it  is  difficult  to  assign 
physical  meaning  to  its  individual  components  (i.e.  the  Dij).  However,  when  the  damage  tensor 
is  expressed  in  its  principal  coordinate  system,  the  eigenvalues  of  D  (denoted  Di,  D2,  and  D3) 
have  a  simple  physical  interpretation.  The  i*  eigenvalue  D,  represents  the  fractional  reduction 
in  load  carrying  area  on  planes  that  are  perpendicular  to  the  i^**  principal  direction.  Thus  the 
use  of  a  symmetric  2"'*  order  damage  tensor  implies  that  the  microscopic  damage  is  orthotropic 
(i.e.  a  system  of  microcracks  and  microvoids  that  are  oriented  along  three  mutually 
perpendicular  directions.  This  type  of  damage  tensor  is  most  often  appropriate  for  fiber- 
reinforced  composite  laminae  since  the  predominate  modes  of  microscopic  damage  are  fiber 
breakage,  fiber  buckling,  fiber/matrix  debonding,  and  matrix  cracks  that  are  oriented  either 
parallel  to  the  reinforcing  fibers  or  perpendicular  to  the  reinforcing  fibers.  Furthermore,  this 
orderly  arrangement  of  microcracks  can  be  described  to  an  acceptable  level  of  accuracy  by 
using  a  2"**  order  damage  tensor  D  whose  principal  directions  are  assumed  to  coincide  with  the 
principal  material  directions  of  the  lamina  (Barbero  and  DeVivo  [11]). 

Damage  tensor  expressed  in  Damage  tensor  expressed  in 

Arbitrary  Coordinate  System  Principal  Material  Coordinate  System 


■  Dll 

Di2 

Di3 

■  Dn 

0 

0  ■ 

■  D, 

0 

0  ■ 

D2I 

D22 

D23 

-> 

0 

D22 

0 

or 

0 

D2 

0 

_  D3I 

D32 

D33  - 

.  0 

0 

D33  - 

.  0 

0 

D3  . 

This  is  a  particvilarly  convenient  assumption  since  it  permits  us  to  work  directly  with  the 
physically  meaningful  eigenvalues  of  D  in  a  coordinate  system  that  we  already  know  from  the 
material  description.  The  eigenvalues  of  the  damage  tensor  must  be  in  the  range  0  ^  Di  <  1, 
where  Di  =  0  corresponds  to  the  case  of  complete  lack  of  microcracks  that  are  normal  to  the  i* 
principal  material  direction,  while  Di  =  1  corresponds  to  the  case  of  complete  separation  of  the 
material  across  planes  that  are  normal  to  the  i*  principal  material  direction.  In  practical 
terms,  the  material  becomes  unstable  and  ruptures  when  the  damage  reaches  a  critical  value 
that  is  most  often  considerably  less  than  1.  Thus  the  practical  ranges  for  the  three  eigenvalues 
of  the  damage  tensor  are  expressed  as 

0  <  Di  <  Dr  i  =  1,2,3  ,  (2) 

cirit 

where  the  values  Di  <  1  (i  =  1,2,3)  are  considered  to  be  material  properties  indicative  of  the 
fracture  toughness  of  the  material  in  planes  normal  to  the  principal  material  axes. 

Related  to  the  damage  tensor  is  the  inteprih/ tensor  £2  =  DgeiCj  given  by 

ft-n  sI-D,  or  Q  =  Vl-D  Qi  =  >/  1  -  Di  i=  1,2,3,  (3) 

where  Qi  and  Di  are  the  eigenvalues  of  the  integrity  tensor  and  the  damage  tensor 
respectively.  Thus,  while  the  eigenvalues  of  the  damage  tensor  D  provide  a  measure  of  the 
fractional  reduction  in  the  original  load  carrjdng  area,  the  eigenvalues  of  the  integrity  tensor  £2 
provide  a  measure  of  the  fraction  of  the  original  area  that  is  still  available  to  cany  load.  Due 
to  the  definition  of  the  integrity  tensor  in  Eq.  (3),  the  principal  directions  of  the  integrity  tensor 
also  coincide  with  the  principal  material  directions,  and  the  eigenvalues  of  the  integrity  tensor 
£2i,  £22,  and  £23  represent  the  square  roots  of  the  available  fraction  of  the  original  load  canying 
area  on  planes  that  are  perpendiciilar  to  the  1,2,  and  3  directions  respectively. 
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1.6  Damaged  constitutive  relations  (stiffness  reduction  scheme) 

The  damage  tensor  D  and  the  integrity  tensor  SI  can  be  used  to  define  the  concept  of  effective 

stress  a  .  For  damaged  materials,  the  definition  of  the  effective  stress  a  represents  an  attempt 
to  assign  the  actual  internal  material  forces  to  cross  sectional  areas  that  have  been  reduced 

due  to  the  presence  of  damage.  Qualitatively  speaking,  we  expect  the  effective  stress  a  to  be 
greater  than  the  apparent  stress  a  (which  uses  the  original  undamaged  cross  sectional  area); 
however,  the  precise  manner  of  defining  the  effective  stress  is  somewhat  arbitrary,  and 
consequently  numerous  methods  have  been  proposed  in  the  literature  (Skrzypek  and 
Ganczarski  [10]). 

In  this  study,  we  use  the  symmetric  effective  stress  tensor  described  by  Cordebois  and 
Sidomoff  [121,  which  can  be  expressed  as 

a  =  m'So  =  (ft’*®  ft’*):o  =  (  (Vl  -  d)'*  ®  (Vl  -  d)’*  )  :o  ,  (4) 

where  M  is  the  4^^  order,  doubly  symmetric,  damage  effect  tensor  whose  components  are  Myki  = 
QikQji.  Thus  as  damage  accumulates,  (i.e.  as  the  Di  increase,  or  as  the  Qi  decrease)  the  load 
carrying  area  decreases  which  causes  the  effective  stress  to  increase  above  the  apparent  or 
nominal  stress.  When  the  above  relationship  is  expressed  in  the  principal  material  coordinate 
system,  then  the  2^^  order  damage  tensor  D  and  the  order  integrity  tensor  Cl  can  be 
expressed  in  terms  of  their  eigenvalues: 


^ 

fti 

0 

0 

Q2 

I 

0  0 

,  and  D 

~  Dx 
0 

0 

D2 

1 - 

0  0 

.  0 

0 

fts  J 

.  0 

0 

D3  j 

Further,  if  we  represent  the  effective  and  apparent  stress  tensors  using  contracted  notation, 

then  we  can  represent  the  damage  effects  tensor  M  and  its  M  ^  inverse  in  the  principal  material 
coordinate  system  as  a  diagonal  6x6  matrix  as  follows. 


M  ^ 


ft?  0  0  0  0  0 

0  ft?  0  0  0  0 

0  0  ft?  0  0  0 


0  0 
0  0 
0  0 


0 

0 

0 


ftgfta 

2 

0 

0 


(5) 
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0 


0 


0 


M 


-1 


0 

0 

0 

0 


0 


^2 


0 

0 


0 

0 

0 

2 

QalQi 


0  0  0  0 


0 

0 

0 

0 

2 

Q1Q2  — 


(6) 


This  jdelds  the  following  expressions  for  the  components  of  the  effective  stress  tensor  o 
expressed  in  the  principal  material  coordinate  system. 


04 

^6 


Q? 


£22 


£2^ 


£22£23 


£23£2i 

_£§_ 

£2i£22 


(l-Di)  ’ 

^2 

(1-D2) ' 

(I-D3)  ’ 

Vi  -D2V1 -D3  ’ 

Os 

Vl  -D3V1  -Di  ’ 


(7a) 

{7b) 

(7c) 

(7d) 

(7e) 

(7f) 


We  can  now  define  the  components  of  the  effective  strain  tensor  e  by  insisting  that  the  actual 
damaged  configuration  and  the  fictitious  undamaged  effective  configuration  should  exhibit  the 
seime  form  for  the  strain  energy  density.  This  is  the  so-called  principle  of  equivalent  strain 
energy,  which  can  be  expressed  as  follows 


1  -  -  1 

2a:e  = 


Now,  substitute  the  inverse  of  Eq.  (4),  i.e.,  a  =  M  :  a  . 

a:e  =  (M:a):e  e:a=e:M:a.  (9) 

Equation  (9)  must  be  true  for  arbitrary  a,  which  yields 


E  =  e :  M  ~  M  ;  e  .  (10a) 

Thus,  based  on  the  principle  of  equivalent  strain  energy,  the  components  of  the  effective  strain 
tensor  e  are  given  by 
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e  =  M  ;  e  =  (Q  0  d)  0(^/1 -  d)  :  e  ,  (10b) 

where  Myki  =  When  expressed  in  the  principal  material  coordinate  system,  the  tensors  D, 

Q,  and  M  can  be  represented  by  diagonal  matrices,  and  thus  yield  the  following  expressions  for 
the  strain  components. 


EjQ?  =  ei(l  -  Di)  , 

(lla) 

CjQa  =  62(1-02)  , 

(11b) 

£3^23  =  £3(1  —  D3)  , 

(11c) 

Z^Cl2^3  ~  £4^  1  —  D2  V 1  “  1^3  > 

(lid) 

egQsOi)  =  egVl  -  D3 'sjl  -Di , 

(lie) 

e6£2iQ2)  =  EgV  1  -  Di  V 1  -  Da  . 

(Ilf) 

The  presence  of  distributed  damage,  as  defined  by  the  internal  variable  D,  affects  the 
Helmholtz  free  energy  'F,  which  is  postulated  to  be  the  sum  of  elastic  strain  energy  (p  and 
dissipation  energy  n: 

¥  =  ^(e,  D,  5)  =  9(e,  D)  +  n(5)  .  (12) 

In  Eq.  (12),  we  have  assumed  introduced  another  internal  variable  5,  which  is  a  non- 
dimensional  parameter  that  will  be  shown  later  to  govern  the  evolution  of  the  damage  tensor  D. 

Assuming  that  the  homogenized  composite  material  exhibits  a  linear  elastic  stress /strain 
relationship  up  to  the  point  of  damage  initiation,  we  can  express  the  strain  energy  q)  in  either  a 
strain-based  form  or  stress-based  form  as  follows 


q)(e,  D)  =  e:C:e  =  ^  are 

(p(a,  D)  =  2  ^  ^  ^  "  2 


(13a) 

{13b) 


where  C  =  C(D)  is  the  4^  order  damaged  material  elasticity  tensor,  and  C  ^  =  C  ^D)  is  the  4*^ 
order  damaged  material  compliance  tensor.  The  dependence  of  these  two  tensors  on  the 
damage  tensor  D  is  defined  by  requiring  that  the  elastic  strain  energy  of  the  actual  damaged 
configuration  should  be  equivalent  to  the  stain  energy  of  the  fictitious,  undamaged  effective 
configuration.  This  strain  energy  equivalence  can  be  expressed  as  follows: 

cp(a,D)  =  cp(a,D=0)  stress-based  (14a) 


or 

(p(e,D)  =  (p(e,D=0)  strain-based  (14b) 


Equations  (14a)  and  {14b)  can  be  written  as 


2^  ‘C  :a  =  ^a:C:a, 

I  I  ^  ^  ^ 

2 €:C: e  =  2 ^  ^  ^ 


(15a) 

(15b) 


where  C  and  C 


represent  the  stiffness  and  compliance  tensors  of  the  undamaged  material. 
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Substituting  Eq.  (4)  a  =  M  ^  0  into  Eq.  (15a)  and  substituting  Eq.  (10)  e  =  M;  e  into  Eq. 

(15b),  and  further  noting  that  a  :  M  ^  M  \*0  and  e  :  M  =  M  :  e  since  0  and  e  are  both 
S3nnmetric  tensors,  we  arrive  at 


2  <T  :  C  :  <j 

=  M M’':a  , 

(16a) 

1  ^ 

2  e:C:  e  = 

^  e  :  M  :  C:  M  :  8  . 

(16b) 

Equations  (16a,b)  must  be  valid  for  any  arbitrary  0  and  e.  Therefore,  Eq.  (16a)  3delds  the 
following  expressions  for  the  damaged  compliance  tensor 

c"^  =  , 

c"'  =  :C'*:  ft"')  , 

c"'  =  ( 0  )  tc"':  ( 0  ) , 

and  Eq.  (16b)  yields  the  following  expressions  for  the  damaged  stiffness  tensor 


M:C:M  ,■ 

(18a) 

(ft  0  ft)  :  C  :  (ft  0  ft)  , 

(18b) 

(  (Vl  -  d)  0  (Vl  -  d)  )  :  C  :  (  (>/l  -  d)  0  (^/l  -  d)  ). 

(18c) 

Eqs.  (17)  and  (18)  define  the  damaged  4^  order  compliance  tensor  C  ^and  the  damaged  4^ 
order  elasticity  tensor  C  in  terms  of  the  4^  order  damage  effects  tensor  M,  the  2^^  order 
integrity  tensor  O  and  the  2^^  order  damage  tensor  D. 

Expressing  Eqs.  (18a-c)  in  the  principal  material  coordinate  system  and  using  contracted 
notation,  the  resulting  4^  order,  damaged  material  elasticity  tensor  C  can  be  represented  by 
the  following  6x6  matrix  which  is  indicative  of  an  orthotropically  damaged  material. 

Cii  C12  Ci3  0  0  0 

C21  C22  C23  0  0  0 

C31  C32  C33  0  0  0 

0  0  0  C44  0  0 

0  0  0  0  Css  0 

0  0  0  0  0  Cee 


The  individual  components  of  the  damaged  orthotropic  material  elasticity  tensor  are  defined  as 


Cii 

=  Ciiftl  =  Cu(l-Di)2  , 

(20a) 

Cl2 

=  Ci2ft?ft2  =  Ci2(1-Di)(1-D2)  =  C21, 

(20b) 

Ci3 

=  Ci3ft?ft3  =  Ci3(1-D,)(1-D3)  =  C31, 

(20c) 

C22 

=  C22ft2  =  C22(l— 02)^  , 

(20d) 

C23 

=  C23ft2ft3  =  C23(1-D2)(1-D3)  =  C32  , 

(20e) 


(17a) 

(17b) 

(17c) 
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C33  =  C33^^3  - 

£33(1-03)2 , 

(20f) 

C44  ==  *2  C44f^2fi3 

=  1  C44(1-D2)(1-D3)  , 

(20g) 

C55  =  ^2  CssOsQi 

=  1  C55(1-D3)(1-Di)  , 

(20h) 

C66  ~  2 

=  ^  C66(1-Di)(1-D2)  . 

(20i) 

Similarly,  the  individual  components  of  the  damaged  orthotropic  material  compliance  tensor 
are  expressed  in  the  principal  material  coordinate  system  as 


Cn 

Cil 

(21a) 

{1-Di)2  ’ 

Cis 

£12 

£12 

C21 , 

(21b) 

(1-Di)(1-D2) 

^'1 

C?3 

£13 

C31  , 

(21c) 

Ci3 

(1-Di)(1-D3) 

C22 

£22  _ 

£22 

(2  Id) 

4 

Q2 

(1-02)2  ’ 

C23 

£23 

£23 

C32  , 

(21e) 

2  2 

Q2O3 

(l-02)(l-03) 

C33 

£33  _ 

£33 

(2  If) 

at 

(1-03)2 

C44 

2C^ 

2CZ 

(21g) 

2  2^ 
Q2  Q3 

(l-02)(l-03)  ’ 

C55 

2655 

2£ss 

(21h) 

2  2 

Q3O1 

(l-03)(l-0j)  ’ 

Cee 

2£66 

2£66 

(21i) 

2  2 

QIQ2 

(1-Oi)(1-D2)  • 

1.7  Damaged  surface 

For  the  independent  kinematic  state  variables  (e,  D,  5)  that  are  used  to  express  the  Helmholtz 
free  energy,  the  associated  energy  conjugate  variables  (i.e.  thermodynamic  forces)  are  obtained 
by  partiad  differentiation  of 'F  with  respect  to  the  kinematic  variables. 

Pairtial  differentiation  of  p'F  with  respect  to  the  strain  tensor  e  yields  the  stress  tensor  a  that 
acts  on  the  damaged  configuration. 

o  =  pa'P/ae  =  C  :  e  =  (m:C:m)  :  e  .  (22) 
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Partial  differentiation  of  T  with  respect  to  the  damage  tensor  D  3delds  the  energy  release  rate 
tensor  Y,  so  named  because  it  represents  the  energy  released  per  unit  increase  of  the  damage 
tensor. 


_  3y(e,D,S)  _  acp(e,D) 

dD  ~  dD 


(23) 


In  Eq.  (24),  C  =  C(D),  while  e  is  independent  of  D. 


Y 

Y 


1  dC  1  d{  M:C:m) 

2  ®  •  3D  •  ®  2  ®  •  aD  •  ®  ’ 

1  ^“dM  3M  „ 

-  2e:M:C:-^:E  =  -e:-^:C:M:e  . 


(24a) 

(24b) 


When  Eq.  (24)  is  evaluated  in  the  principal  material  coordinate  system,  the  energy  release  rate 
tensor  Y  is  diagonal  and  its  eigenvalues  Yi  (i=l,2,3)  are  expressed  in  terms  of  the  actual  strains, 
the  undamaged  stiffness  and  the  eignevalues  of  the  integrity  tensor  as 

Yi  —  CiifiiEi  +  Ci2f^2  £ie2  Ciafita  £163  +  £5  +  C66f^2  £6  ,  (25a) 

Y2  =  C21Q1  £i£2  +  022^^2  £2  C23f^3  £2£3  +  044^^3  Es  Ceef^i  £6  >  (25b) 

Ya  =  CiiQi  E183  +  C32Q2  £2£3  +  C23f^3  £3  +  044^^2  £5  +  C55Q1  £5  .  (25c) 

Substituting  Eqs.  (20a-i)  into  Eqs.  (25a-c),  the  eigenvalues  of  the  energy  release  rate  tensor  Yi 
(i=  1,2,3)  can  be  expressed  in  terms  of  the  actual  damaged  stiffness  tensor,  the  actual  strain 
tensor  and  the  integrity  tensor: 


Yi  = 

1  (  2 

—  ICnei  +  Ci2£i£2  +  Ci3£i£3 

+  C55  £5 

+  CeeEe), 

(26a) 

Ya  = 

1  /  2 

~2  1  ^21  £162  +  C22  £2  +  C23  £2£3 
na 

+  C44  £5 

+  Cee  £6  ) , 

(26b) 

Ys  = 

~2  (  ^11  ^1^3  ^32  £2£3  +  C23  £3 

Q3 

+  C44  £5 

+  C55  £5  ) . 

(26c) 

Note  that  the  potential  (p  in  Eq.  (23)  can  alternately  be  expressed  using  e(D)  and  C,  instead  of 
e  and  C,  which  would  yield 


Y  = 


Y  =  - 


(f 


e  :  C  :  e 


=  -  p3(p/3D  =  - 


dD 


dD 


lde-~ 
2  dD  •  ^  ® 


1  ~  ~  de 

2  ®  ^ • dD 


de  ~  ~ 
dD  •  ® 


(27a) 

(27b) 


When  evaluated  in  the  principal  material  coordinate  system,  Eq.  (27b)  yields  the  same 
expressions  for  the  eigenvalues  of  Y  as  obtained  earlier  in  Eqs.  (25)  and  (26). 

We  can  express  the  energy  release  rate  tensor  Y  in  terms  of  stress  components  by  utilizing  the 
dual  of  the  potential  (p  as  follows 
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o :  e  =  (p(e,D)  +  (p*(CT,D)  ->  (p(e,D)  =  a :  e  -  9*(a,D)  .  (28) 

Thus  the  energy  release  rate  tensor  can  be  expressed  in  either  of  the  following  ways. 

^  “  -PaD  =  '  (29) 


where  (|)*(o,  D)  =  ^  a  :  C  S  a  and  C  *  =  C  *(D). 


Substituting  Eq.  (30)  into  Eq.  (29)  yields 


Y  = 


a(^a;C~Sq) 


ac 


-1 


aD 


Y  - 


n 


"1  —“1 


2  ®  •  aD 
1 

:c  -  20 


:  c  = 


=  C  (D). 

(30) 

1  ^ 

2'^- 

a5 - ^  = 

(31a) 

■1  aM'' 

=  aD 

=  -q:  3D  :M  :q. 

(31b) 

When  Eq.  (31b)  is  evaluated  in  the  principal  material  coordinate  system,  the  eigenvalues  of  the 
energy  release  rate  tensor  Yi  (i=l,2,3)  can  be  expressed  in  terms  of  the  undamaged  compliance 
tensor,  the  actual  stress  tensor,  and  the  integrity  tensor  as  follows. 


Y, 

Yz 

Ys 


fcil  2 

QtU;'” 
fei 


2^ 

^4 


2 


Cl2 


2  <^102  + 


-^0102 

C22  2 
2  ^2 
Q2 


C32 

+  — r  ( 


2  CTi03 

Q3 

C23 
2  0203 
Q3 


C55  2 

2  ^5 
Qi 

+ 

Cee  2 

2  <^6 
02  y 

—1  ^ 

Cee  2 

2  ^6 
Oi  ; 

C44  2 

2 

Qs 

+ 

C44  2 

2 

O12 

+ 

C55  2I 

— 2  ^5 

£2?  J 

(32a) 

(32b) 

(32c) 


Substituting  Eqs.  (21)  into  Eqs.  (31),  the  eigenvalues  of  the  energy  release  rate  tensor  can  be 
expressed  in  terms  of  the  actual  damaged  compliance  tensor,  the  actual  stress  tensor  and  the 
integrity  tensor. 


Yi  “  2  (^11  ^12  O1G2  +  Ci3  O1C3 

+  CsUs 

+  Cee  06  ), 

(33a) 

2  2  -1 
*2  “  2  1^21  CT2CT1  +  C22  02  +  C23  G2O3 

O2 

-1  2 
+  C44  O4 

C66  Os  )  > 

(33b) 

Y3  =  -  ^C31  O3CT1  +  C32  CT302  +  C33  O3 

O3 

+  C44  CT4 

+  Cssas). 

(33c) 

Equation.  (32)  and  (33)  can  also  be  obtained  by  expressing  the  dual  potential  9*  in  terms  of  the 
effective  stress  a  and  the  undamaged  compliance  C  \  i.e.,  9*(0,  D)  =  C  S  q  where 

q  =  q(D)  and  C  ^  is  independent  of  D. 
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The  damage  surface,  which  is  assumed  to  separate  non-damaging  behavior  from  damage- 
inducing  behavior,  is  defined  as  a  complete  quadratic  function  of  the  energy  release  rate 
tensor. 

g(Y,7)  =  ^[Y^  +  VIH:Y|  -  {#)  +  Yo).  (34) 

In  Eqs.  (33),  the  4^^  order  tensor  J  and  the  2^^^  order  tensor  H  represent  experimentally  derived 
material  constants  that  define  the  damage  tolerance  of  the  material.  When  expressed  in  the 
principal  material  coordinate  system,  J  and  H  are  each  assumed  to  have  only  three  non-zero 
values,  denoted  Jn=  Jim,  J22  =  J2222,  J33  ^  J3333,  Hi  =  Hn,  H2  =  H22,  H3  =  H33  .  These  six 
material  constants  define  the  shape  of  the  damage  surface.  The  material  parameter  yo  is  also 
determined  from  experimental  data  and  defines  the  initial  size  of  the  damage  surface,  thus 
functioning  as  an  initial  damage  threshold.  In  Eq.  (33),  the  linear  term  H:Y  is  necessary  to 
define  a  material  that  has  different  damage  tolertoce  in  tension  and  compression.  When 
expressed  in  the  principal  material  coordinate  system,  the  damage  surface  simplifies  to  the 
following  form. 


g(Y,Y)  =  VJuY?  +J22Y^  +J33Y^  +  ViHiYi  +  H2Y2  +  H3Y3I  -  (  Y(8)  +  Yo )  (35) 

Depending  upon  whether  the  energy  release  rate  eigenvalues  Yi  (i=  1,2,3)  are  expressed  using 
Eqs.  (25),  (26),  (32),  or  (33),  the  damage  surface  can  be  defined  in  terms  of  either  of  the 

following  sets  of  variables:  (C,i2,8,')^,  (C,^2,e,')^,  (C  \i2,a,y),  or  (  C  \Q,a,y)  . 

The  damage  surface  defined  by  Eq.  (35)  has  the  shape  of  the  Tsai-Wu  failure  surface  when  it  is 
expressed  in  stress  space.  This  can  be  seen  by  noting  that  the  energy  release  rate  eigenvalues 
[Eqs.  (32)  or  (33)]  are  quadratic  functions  of  the  stress  components,  thus  the  damage  surface 
[Eq.  (35)1  can  be  expressed  as 

g(Y,-))  =  PyCiOj  +  ttiOi  -  (Y  +  Yo)  i,j  =  1,2,  -,6  ,  (36a) 

where  the  coefficients  Py,  and  ai  are  functions  of  the  damaged  compliance  tensor  C~i(D),  the 
integrity  tensor  12(D),  and  the  material  damage  characteristic  tensors  J  and  H.  The  Tsai-Wu 
failure  criterion  [13]  predicts  failure  when 

FijOiOj  +  FiOi  >1  i,j  =  1,2,. ..,6  ,  (36b) 

where  the  coefficients  Fy  and  Fi  are  functions  of  the  lamina  strength  in  the  various  modes  of 
deformation. 

If  we  choose  to  define  the  quantities  y  eind  yb  such  that  the  expression  (  y  +  yo  )  =  1  indicates 
complete  failure  of  a  lamina,  then  the  damage  surface  g(Y,'y^  increases  failure  in  the  same 
manner  as  the  Tsai-Wu  criterion  [13]. 

The  damage  surface  [Eq.  (35)]  defines  the  limits  of  non-damaging  behavior  at  an  arbitrary  point 
in  the  lamina.  If  the  current  local  values  of  Yi  (i=l,2,3)  and  y  yield  a  point  on  the  surface 
g(Y,y),  then  further  damage  is  assumed  to  occur.  The  damage-flow  surface,  denoted  as  f(Y,'^, 
is  defined  similar  to  the  damage  surface  g(Y,y)  but  without  the  linear  term  H:Y. 

f(Y,Y)  =  V  Y:  J:Y  -  (y(5)  +  Yo  )  =  -  (  Y(5)  +  Yo)  •  (36c) 
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1.8  Damage  evolution 

The  dissipation  potential  n  is  assumed  to  be  of  the  form 
jc(8)  =  -ci(^8  +  C2exp(^^^^  , 


(37) 


where  the  dimensionless  parameter  8>0  controls  the  evolution  of  the  damage  process.  This  is 
the  simplest  dissipation  potential  with  the  minimum  number  of  adjustable  coefficients  that  still 
captures  the  experimentally  observed  behavior  of  fiber-reinforced  composite  laminae.  The 
evolution  law  is  assumed  to  be  isotropic  for  simplicity  and  due  to  the  lack  of  experimental 
observations  indicating  anisotropic  evolution  of  the  damage  surface  (Barbero  and  DeVivo  [11]). 

The  energy  conjugate  variable  to  5  is  denoted  as  y  and  is  obtained  by  partial  differentiation  of 
the  free  energy  with  respect  to  5: 


(38a) 


The  evolution  of  the  damage  tensor  D  occurs  in  a  direction  that  is  normal  to  the  damage-flow 
surface  as  defined  by  Eq.  {36c): 


J:Y 


+J22Y2  +J33’^ 


JiiYi  +  J22Y2  +  J33Y3 


(38b) 


The  evolution  of  the  damage  surface  is  given  by 


g  =  =  p(-l)  =  -p.  (38c) 

In  the  interest  of  keeping  the  size  of  the  report  in  reasonable  limit,  numerical  results  based  on 
the  damage  models  discussed  herein  are  not  reported.  They  will  appear  in  the  literature  in  the 
near  future. 


2.  Finite  Element  Formulation  of  the  Layerwise  Theory 
2.1  Introductory  comments 

The  layerwise  theory  proposed  by  Reddy  [14-16]  exploits  the  laminate  architecture  to  separate 
the  thickness  variation  from  the  surface  variations  and  the  thickness  variation  is  represented 
independently  using  Lagrange  polynomials.  Reddy  and  Mitchell  [17]  presented  a  general 
formulation  of  the  problem  using  both  ESL  and  layerwise  theories.  Layerwise  finite  element 
models  of  Reddy's  layerwise  theory  can  be  found  in  [14-16].  The  work  of  Saraveinos  et  al.  [18]  is 
perhaps  the  first  one  to  use  a  layerwise  finite  element  model  for  the  analysis  of  piezoelectric 
plates.  However,  these  authors  only  considered  constant  aind  linear  approximation  schemes 
across  the  thickness. 

In  the  present  study,  a  general,  displacement  finite  element  formulation  of  the  layerwise  theory 
of  Reddy  [14-16]  for  laminated  plate  structures  with  piezoelectric  materials  (layers  or  patches) 
is  developed.  The  formulation  includes  full  electromechanical  coupling  and  arbitrary 
approximation  through  the  thickness  as  well  as  in  the  surface  of  the  laminate.  Several 
approximations  are  used  for  the  primary  variables  of  the  problem  in  the  thickness  direction 
and  different  interpolation  schemes  are  considered  in  the  surface  directions. 
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2.2  Theoretical  formulation 

Consider  a  laminated  plate  with  thickness  H  and  arbitrary  planar  geometiy,  built  with 
piezoelectric  laminae  or  patches  and  laminae  made  with  different  linear  elastic  materials.  A 
global  orthonormal  recteingular  reference  frame  (x,i/,z)  with  the  z  axis  aligned  with  the  laminate 
thickness  is  assumed.  The  material  within  each  layer  of  the  laminate  is  assumed  to  be 
homogeneous,  generally  anisotropic  and  lineair  elastic,  including  linear  piezoelectric  effects.  The 
constitutive  law  considered  for  a  composite  materiad  is  [19-21]: 

Di=xfhi  (39) 


and  for  a  piezoelectric  material: 


where  Sij  and  eim  are  the  components  of  the  stress  and  strain  tensors,  Di  and  Ei  are  the  electric 
displacement  and  electric  field  vector  components,  Cijim  are  the  components  of  the  fourth-order 
tensor  of  elastic  moduli,  euj  are  the  components  of  piezoelectric  moduli  and  Xii  are  the 
components  of  the  dielectric  moduli  for  the  /c^  lamina.  The  constitutive  behaviour  of 
piezoelectric  materials  makes  a  piezolaminated  structure  a  coupled  mechanical  and  electrical 
problem. 


The  equilibrium  of  stresses  at  the  interface  between  adjacent  plies  (/c)  and  {k  + 1)  requires  the 
stress  field  to  satisfy  the  conditions  (see  Figure  1) 


(^+1) 


(41) 


The  electric  displacements  at  the  interface  must  satisfy: 

/)(*-•)  .  )  =  w. 


(42) 


where  We  is  the  imposed  surface  charge.  If  no  electric  potential  is  prescribed  along  the 
interface,  or  equivalently  no  surface  charge  is  imposed,  the  electric  displacement  Dz  component 
will  be  continuous  across  the  interface.  Because,  in  general,  the  material  properties  differ  for 
layers  {k)  and  (A:  +  l) ,  conditions  (41)  and  (42)  can  only  be  satisfied  if: 


(^+1) 


(43) 


These  conditions  imply  that  the  displacement  and  electric  fields  must  be  continuous  in  the 
thickness  direction,  ^ong  the  interfaces  joining  dissimilar  materials.  Hence,  the  primary 
variables,  neimely,  the  displacement  field  components  (u,  w)  and  electric  potential  cp,  are 
represented  as  products  of  functions  of  and  C®  continuous  functions  of  the  thickness 

coordinate  z  (the  upper  indices  are  not  exponents): 
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smart  layer 


Figure  1:  Equilibrium  of  interlaminar  stresses. 


^  I  T 

u(x,y,z,t)='^u\x,y,t)F^ (z);  v{x, y,z,t)='^^  (^> y> 0^  (z) 

'Si  't 

wix, y,z,t)='^W^ (x, y, t)Y^ (z);  jix, y,z,t)='YjJ^ Y’ (^) 
7=1  7=1 


Following  the  approach  proposed  by  Reddy  [14],  the  thickness  direction  continuous 
functions  Y^(2!j  and  (y(z)  are  constructed  using  the  Lagrange  pol3niomials  and  have  the 

following  characteristics: 


/  T  I 

(a)  P\^j)  =  dij.  Y\zj)  =  djj.  Q\zj)  =  dij.  where  ^IJ  =\  (45) 

[O,  7^7 

(b)  local  support,  the  functions  are  nonzero  only  in  the  layers  that  share  the  Rh  interface 
(see  Figure  2). 

With  this  approach  several  finite  element  models  can  be  constructed  in  an  efficient  and 
straightforward  manner. 

The  unknowns  for  the  this  layerwise  theory  are  the  surface  functions  V^(jc,2/,1),  W^(x,y,f) 

and  j^{x,yjtj.  The  number  of  functions  Nf  ,  Ny  and  Nq  ,  considered  by  the  theory,  depend  on  the 
number  of  layers  and  on  the  degree  of  the  assumed  approximation  along  the  thickness  of  each 
layer  for  each  one  of  the  primsoy  variables.  Underlying  this  assumption  for  the  primary 
variables  is  the  hypothesis  that  the  behaviour  of  these  fields  along  the  surface  directions  (x,z/), 
is  different  from  the  behavior  along  the  thickness  direction  (2^.  This  is  an  assumption  common 
to  all  the  plate  and  shell  theories,  that  allows  reducing  a  three-dimensional  problem  to  a  two- 
dimensional  one. 
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Figure  2.  (a)  Layer-wise  displacement  field,  (b)  Laywerwise  representation  of  a  laminate. 


The  layerwise  finite  element  models  are  developed  using  a  variational  formulation.  For 
details  on  the  theory  and  veiriational  formulation  of  the  layerwise  theory,  the  reader  may 
consult  References  14  through  16.  To  obtain  the  finite  element  model,  the  primary  variables 
are  approximated  inside  each  element  as: 


The  numbers  W,  N}/,  Nw^  and  N/  of  interpolation  functions  aK^(x,y,tj,  CK^(x,yytl  and 

depend  of  the  interpolation  scheme  associated  with  the  Ah  layer  thickness  function 
inside  the  element.This  set  of  equations  can  be  cast  in  a  matrix  form: 
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After  assembling  the  contributions  from  all  the  elements,  the  global  system  of  equations 
obtained  can  be  written  as: 

[M]{A}+[/i:]{A}  =  {F},  (48) 

where  [Af|  is  the  mass  matrix,  [K]  the  stiffness  matrix,  (J^)  the  load  vector,  and  {D}is  the  vector  of 
nodal  values  of  the  generalized  displacements. 

Numerical  results  of  the  formulation  presented  above  can  be  found  in  [23]. 


19 


3.  Finite  Element  Analysis  of  Plates  with  Embedded  ActuatorsUsing  FSDT 

3.1  Introduction 

The  displacement  fields  of  all  equivalent  single-layer  plate  theories  are  expanded  in  terms  of  the 
thickness  coordinate  up  to  a  desired  degree.  However,  theories  higher  than  the  third  are  not 
attempted  due  to  very  little  improvement  in  accuracy  for  the  amount  of  algebraic  complexity 
and  computational  effort  involved.  Here  we  discuss  the  nonlinear  finite  element  formulation  of 
the  first-order  shear  deformation  plate  theory  (FSDT),  which  will  be  used  in  Phase  II  effort  in 
conjunction  with  the  layerwise  theory  of  Reddy  [14,15]  to  carry  out  the  global-local  analysis  of 
composite  laminates. 

3.2  Equations  of  motion 

The  displacement  field  for  the  first-order  shear  deformation  theory  can  be  expressed  in  the  form 
(see  Reddy  [15]) 

«ix,y,z,t)  =  Uo(x,y,t)  + 

v(x,y,z,t)  =  Voix,y,t)  +  z(py(x,y,t)  (49) 

w{x,y,z,t)  =  Wo{x,y,t) 


where  (wo>^o>^o)  <ienote  the  displacements  and  rotations  of  transverse  normals  on 


the  plane  z  =  0,  respectively  (see  Figure  3). 


Figure  3.  The  kinematics  of  deformation  of  a  typical  edge  in  the  first-order  shear 
deformation  plate  theory  (FSDT). 

The  equations  of  motion  of  the  first-order  shear  deformation  theory  are  derived  using  the 
dynamic  version  of  the  principle  of  virtual  displacements  (see  Reddy  [15]).  The  equations  of 
motion  are 
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are  the  mass  inertias,  and  (N^.NyyyN^)  denote  the  total  in~plane  force 
(M^.MyyyM^ )  the  moment  resultants: 
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The  force  and  moment  resultants  are  related  to  the  strains  by 
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Here  e^j  are  the  transformed  moduli  of  the  actuating/ sensing  material,  and  is  the  electric 
field.  The  coefficients  of  A^j ,  5,yand£)^  are  given  in  terms  of  the  layer  transformed  stiffnesses 
and  layer  coordinates  and  z^  (see  Figure  4). 


Figure  4.  Schematic  of  a  laminated  plate. 


3«3  Velocity  feedback  control 

Here  we  consider  the  velocity  proportional  closed-loop  feedback  control,  the  magnetic  field 
intensity  H  is  expressed  in  terms  of  coil  current  I{x,yj)  as 

H{x,y,t)  =  k^I{x,y,f).  (56) 

The  current  lix^yA)  is  related  to  the  transverse  velocity  of  the  plate  wq  by 

Iix,y,t)  =  c{t)^.  (57) 

ot 

Here  denotes  the  coil  constant,  which  can  be  expressed  in  terms  of  the  coil  width  ,  coil 
radius  and  number  of  turns  in  the  coil  by 


(58) 


and  c{t)  is  the  control  gain.  In  this  study  the  control  gain  is  assumed  to  be  constant.  Hence  it 
is  represented  as  c. 
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3.4  Weak  formulation  of  the  equations 

The  weak  forms  of  the  equations  of  motion  are 


d&o  dduQ 

+<S<0 


3^0  Ttr  -KT  c 


J^\ 


bcdy  -  (|^  {dto  +  N^Tjy  )}ds  , 
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+  SwolQ^-^‘dKify-<i^QV„ds, 

a/2  ^r 


xy 


-+N 


dx  ^ 


(59) 


where  V„  is  defined  as 


V„  =  iQ^n,^Q^n^)+{N^^  +  N^^n,+\N^^ 


3.5  Finite  element  model 


^  dy 


\ny. 


(60) 


The  primeuy  variables  of  the  first-order  theory  are  i/^,  Wo,  <f>^,  where  {u„,uj  denotes  in¬ 
plane  normal  and  tangential  displacements,  and  are  the  rotations  of  a  transverse  line 

about  the  in-plane  normal  and  tangent.  Lagrange  interpolation  functions  of  (u„,u^,Wo, 
are  used  for  the  formulation  of  the  displacement  finite  element  model.  The  finite  element 
model  has  five  degrees  of  freedom  per  node  and  four-node,  eight-node,  or  nine-node 
rectangular  elements  may  be  used. 

The  generalized  displacements  £tre  approximated  over  a  typical  element  Q®  by  the  expressions 


m  m 

Uxix,y,t)  =  '^u^^(t)  \|/f(A:,y),  Vx(,x,y,t)  =  '^v^(t)  \|/f(jc,y), 
/=1  /=1 


m  __ 

wo  (x,  0  =  X  > 

/=1 

m  m 

(x,  0 = X  V/  (^.  >') .  <i'y  (x,  T,  0 = X  'I'f  y'^  > 

i=l  /=1 


(61) 
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where  y/^  and  (pf  denote  the  Lagrange  interpolation  functions.  Here  we  chose  the  same 
approximation  for  the  displacements  (Wo3^o>^o)  rotations  although  one  could  use 

different  approximations  for  these  two  sets  of  fields. 

The  finite  element  model  is  of  the  compact  form 


5 

I 

p=i  y=i 


I  *^.7 
•  1 


ap.P 


J  V 


l-fi«=0  ,  1  =  1,2, 


(62) 


where  «  =  1,  2, 3,  4,  5  and  i\=n2  =  n^  =  n^  =  n^.  The  nodal  values  are 

A*,.  =  Uj,  =  Vj,  A)  =  wj.  A)  =  Xj,  A)  =  Yj  (63) 

and  the  non-zero  stiffness,  mass,  and  damping  coefficients  are  defined  in  Appendix  I. 

This  completes  the  general  non-linear  finite  element  model  development  for  the  first-order 
shear  deformation  plate  theory  with  actuator  layers.  The  finite  element  model  is  called  a 
displacement  finite  element  model  because  it  is  based  on  equations  of  motion  expressed  in 
terms  of  the  displacements,  and  the  generalized  displacements  are  the  primary  nodal  degrees 
of  freedom. 

Transient  Analysis 

Here  we  discuss  the  procedures  to  determine  the  transient  response  of  composite  laminates. 
The  equations  of  motion  (62)  must  be  approximated  further  using  the  Newmark  method.  First, 
note  that  Eq.  (62)  has  the  general  form 

[M*]  {«"} +[C^]  m+[K^]  {u^}  =  {F^} 

The  global  displacement  vector  {«}  is  subject  to  the  initial  conditions  that  the  displacement 
and  velocity  are  known  at  time  ^  =  0 :  {u(0)}  -  {«}o  {w(0)}  =  {m}o  .  In  Newmark’s  scheme,  the 

function  (of  time)  and  its  first  derivative  are  approximated  using  Taylor’s  series  and  only  terms 
up  to  the  second  derivative  are  included: 

Ms+1  =M5' +^{“>5 +^(A/)^{«}s+Y’  Ws+1  ={«}5+Mw}5+a.  {w}5+a  =0-a){«}5 +a{«}5+l  (65) 


where  At  is  the  time  increment,  Af  =  -t^  and  is  the  current  time  and  is  the  next 

time  at  which  we  seek  the  solution.  We  assume  that  the  solution  at  time  t^  is  known. 
Substituting  the  third  equation  into  the  first  two  and  solving  for  {m}  ,  we  obtain 

{"h+l  ={“}^  +«l{wh  +^2  {«},+!  .  {«Ul  =a3({“Ul 


r 

and  a  and  y  are  parameters  that  determine  the  stability  and  accuracy  of  the  scheme  Finally 
,we  obtain  the  finite  element  equations  of  the  form 


2  (1 
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={F}s;s.l 


(68) 
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3.6  Numerical  results  -  linear 

Numerical  studies  are  carried  out  to  analyze  the  vibration  suppression  characteristics  for  the 
different  lamination  schemes  and  boundary  conditions  using  the  previously  developed  TSDT 
finite  element  models.  The  baseline  of  the  simulations  is  the  simply  supported  square  laminate 
(a/b  =  l ,  a/h  =  lO)  under  sinusoidal  distributions  of  velocity  initial  conditions  over  the  domain, 

^(x,y,/  =  0)  =  sin— sin— .  The  selected  time  step  in  the  present  study  is  O.OOOSsec.  The 
dt  a  b 

lamination  scheme  (0i,02»03»04,w)s  means  that  the  10  layered  symmetric  laminate  with  the 
fiber  orientation  being  ( 0i, 02, 03, 04, 04,03,02,01),  where  stands  for  the  magnetostrictive 
layer  and  subscript  s  stands  for  symmetric.  The  material  properties  of  Terfenol-D  and  the 
elastic  composite  materials  are  listed  in  Table  1. 

In  finite  element  analysis,  solution  symmetries  should  be  taken  advantage  of  identifying  the 
computational  domain  because,  they  reduce  computational  efforts.  For  a  laminated  composite 
plate  with  all  edges  simply-supported  or  clamped,  a  quadrant  of  the  plate  may  be  used  as  the 
computational  domains.  Figure  5  shows  the  two  types  of  simply  supported  boundary  condition 
for  third-order  shear  deformation  theory.  Quarter  plate  models  with  proper  boundary 
conditions  can  be  used  in  the  antisymmetric  laminates  with  simply  supported  boundary 
condition  but  not  with  the  clamped  cases.  For  symmetric  laminates,  the  simply  supported 
cross-ply  laminates  can  be  modeled  as  a  quarter  plate. 

Table  1.  Material  properties  of  magnetostrictive  eind  elastic  composite  materials. 


Material 

El 

[GPa\ 

E2 

[GPa\ 

Gi3 

\GPa\ 

G23 

\GPa] 

Gi2 

[GPa\ 

P 

[Kgm~^] 

A 

Terfenol-D 

26.5 

26.5 

13.25 

13.25 

13.25 

0 

9250 

1.67 

CFRP 

138.6 

8.27 

4.96 

4.12 

4.96 

0.2 

6 

1824 

- 

Gr-Ep(AS) 

137.9 

8.96 

7.20 

6.21 

7.20 

0.3 

0 

1450 

- 

Gl-Ep 

53.78 

17.93 

8.96 

3,45 

8.96 

0.2 

5 

1900 

- 

Br-Ep 

206.9 

20.69 

6.9 

4.14 

6.9 

0.3 

0 

1950 

- 
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FSDT:  Uo=  1^0=  <t>x=0 


TSDT  &CLPT: 
W0“  U^o  =  W^O.*=0 


Figure  5.  Simply  supported  boundaiy  conditions,  SS-land  SS-2,  for  laminates  using  the  third- 
order  shear  deformation  theory. 


Simply  supported  composite  laminates 

To  compare  with  the  analytical  results,  the  SS-1  boundaiy  conditions  and  quarter  plate  F.E. 
model  have  been  used  in  tiie  symmetric  cross-ply  launmates.  The  vibration  suppression  effects 
from  the  ainal3d:ical  and  numerical  methods  are  within  the  reasonable  range  of  distribution  as 
shown  in  Figure  6.  Figure  7  shows  the  central  displacements  using  the  different  plate  theories 
(CLPT,  FSDT,  and  TSDT)  for  two  different  lamination  schemes.  It  is  observed  that  CLPT  gives 
the  higher  vibration  suppression  capacity  in  both  cases.  This  is  expected  because  the  CLPT 
renders  the  plate  stiffer  compared  to  the  other  theories.  After  studying  the  influence  of  lamina 
material  properties  on  the  amplitude  of  vibration  and  vibration  suppression  times,  it  is 
observed  that  materials  having  the  almost  same  EjE^  ratio  have  similar  vibration 

suppression  characteristics.  Figure  8  shows  the  vibration  suppression  characteristics  for  the 
different  laminate  materials. 


Time  (sec) 


Figure  6.  Comparison  of  analytical  and  numerical  central  displacements  with  the  case  of 
uncontrolled  behaviour  for  the  S3mimetric  cross-ply  CFRP  laminates  (m,90,0,90,0)a  . 
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0.0006 
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Figure  7(a)  The  comparison  of  the  central  displacements  by  the  different  plate  theories  for 
CFRP  simply  supported  laminates;  S3rmmetric  cross-ply  CFRP  laminate  (m,90,0,90,0)s. 


Time  (sec) 

Figure  7(b).  The  comparison  of  the  central  displacements  by  the  different  plate  theories  for 
CFRP  simply  supported  laminates;  symmetric  cross-ply  CFRP  laminate  (0,90,0,90,m)s  . 


Figure  8.  The  effect  of  the  laanina  material  properties  on  the  vibration  suppression 
characteristics  for  the  symmetric  cross-ply  laminates  (m,90,0,90,0)s . 
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Clamped  composite  laminates 

The  same  laminated  plates  but  with  fully  clamped  boundary  conditions  are  analyzed  using  a 
8X8  mesh  in  a  full  plate.  The  effect  of  the  boundary  conditions  on  the  vibration  characteristic 
is  shown  in  the  Figure  9.  The  maximum  displacements  of  the  simply  supported  plate  are 
greater  than  those  of  the  clamped  case,  which  is  expected.  For  the  vibration  suppression,  the 
cases  of  simply  supported  laminates  that  has  the  bigger  displacements  take  less  suppression 
time. 


Time  (sec) 


Figure  9.  The  effect  of  boundary  conditions  on  the  central  displacements  for  the  different 
laminates;  compeirison  of  simply  supported  and  clamped  laminates. 


Effects  of  the  mechanical  loading 


Numerical  studies  are  also  carried  out  to  analyze  the  vibration  suppression  characteristics  for 
the  composite  laminates  under  the  uniformly  distributed  load  instead  of  velocity  initial 

conditions.  Figure  10  shows  the  vibration  suppressions  for  the  selected  simply  supported  and 
clamped  laminated  plates  under  continuously  applied  uniformly  distributed  loading  while 
Figure  1 1  shows  the  case  under  suddenly  applied  and  removed  distributed  loading.  The  effect 
of  sinusoidal  loading  on  the  central  displacement  has  been  studied.  TTie  results  of  the 
symmetric  cross-ply  laminate  with  simply  supported  laminates  subjected  to  the  sinusoidal  and 
uniformly  distributed  loadings  are  shown  in  Figure  12.  Here  the  sinusoidal  loading  is 

sin— sin—  and  the  transverse  displacements  in  Figures  10-12  are  plotted  in  the 
a  b 


nondimensionalized  forms  as  wq  (100) 


E2h^ 
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Figure  10.  Nondimensionalized  central  displacement  for  simply  supported  and  clamped 
laminated  plates  under  uniformly  distributed  loading. 


Figure  11.  Nondimensionalized  central  displacement  for  simply  supported  and  clamped 
laminated  plates  under  suddenly  applied  uniformly  distributed  loading. 


Time  (sec) 


Figure  12.  Comparison  of  central  displacements  for  simply  supported  laminated  plates 
under  Sinusoidal  and  uniformly  distributed  loadings. 
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3.7 Numerical  results  -  nonlinear 

Laminated  composite  square  plate  (alb  =  l)  with  both  the  upper  and  lower  surfaces 
embedded  magnetostrictive  materials,  as  shown  in  Figure  2,  is  considered.  The  plate  is 
made  of  composite  fiber  reinforced  polymer  (CFRP)  composites  and  the  magnetostrictive 
material,  Terfenol-D.  The  adhesive  layers  are  neglected  in  the  analysis.  The  laminated 
composite  plates  are  composed  of  total  10  layers  and  all  the  layers  are  assumed  to  be  of  the 
same  thickness.  Two  side-to-thickness  ratios  a//i  =  10andl00  are  considered  to  represent 
the  thick  and  thin  composite  plates. 

The  geometric  and  material  parameters  used  are 

•  Geometry  :a  =  b  =  l00cm,  h-lOorlcm 

•  CFRP  :  E,  =  138.6  GPa,  =  8.27  GPa,  G^^  =  G,^  =  4.96GPa,  G  23=  4.12GPa 

=  0.26,  p  =  lS24  Kg  I 

.  Terfenol-D  :  E"  =  26.5  GPa,  d”’  =  1.67x1  O'* wJ”',  V  =  0,p  =  9250Kglm^ 

The  numerical  analysis  variables  in  this  study  and  the  representing  notations  that  are 
used  in  the  result  figures  are  the  following: 

•  3  Plate  theories  :  CLPTfQ,  FSOT^Fj,  TSDT(T) 

•  4  Lamination  schemes  :  S5mimetric  cross-ply  (m, 90, 0,90,0)8  fSC),  symmetric 
angle-ply 

(m,45,— 45,45,— 45)s  fSA),  symmetric  general  angle-ply 
(m,45,— 45,0,90)sfSG),  and  asymmetric  general  angle-ply 
(m,45,— 45,15,— 15,0,90,30,  — 30, m)  (AS)  laminated  plates. 

•  4  Loading  conditions  :  Uniformly  distributed  load^UMD),  uniformly  distributed  impact 
load(IMD),  uniformly  distributed  sinusoidal  load(SUMD)f  uniformly  distributed  sinusoidal 
impact  load^SIMD). 

•  The  other  indices:  Linear/L)  and  Nonlinear^JVL)  analyses,  Thick(a/A  =  10)  plateflO)  and 

thm{a/h  =  100)  plate^IOO),  and  Load  intensity  =  10^ (EO),  =  5.0x10^^5^77,  and  so  on. 

The  shear  correction  factor  used  in  FSDT  is  5/6 ;  m  represents  magnetostrictive  layer  and 
subscript  s  stands  for  S3anmetric,  and  the  sinusoidal  load  means 

q{x,  y)  =  q^  smimc/ a)  sminyfE) . 

The  displacements  are  measured  at  the  center  of  the  plates  and  nondimensionalized  as 
^IQOy. E^h^  j a^q^.  A  tolerance  of  5  =  10”^  is  used  for  convergence  in  the  Newton- 
Raphson  iteration  scheme  to  check  for  convergence  of  the  nodal  displacements. 

First,  the  effects  of  the  different  plate  theories  on  the  transient  responses  are  considered. 
Linear  and  nonlinear  responses  obtained  by  the  three  different  theories  are  presented  in 
Figure  13.  It  is  observed  that  the  effect  of  nonlinearity  on  the  tremsient  responses  is  to 
decrease  the  amplitude  and  increase  the  frequency.  Note  that  due  to  the  large  geometric 
nonlinearity  effects  the  nonlinear  transient  behaviors  between  TSDT  and  other  two  theories 
are  apparent.  It  is  also  observed  that  the  CLPT  theory  gives  higher  frequencies  and  lower 
amplitudes.  It  is  because  CLPT  theory  renders  the  plate  stiffer  compared  to  the  other  two 
theories. 
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Next,  the  effect  of  applied  loading  conditions  on  the  vibration  suppression  can  be  seen  from 
the  Figure  14.  The  four  different  loading  conditions  are  considered  to  study  their  effect  on 
the  response.  They  are  uniformly  distributed  load,  uniformly  distributed  sinusoidal  load, 
uniformly  distributed  impact  load  and  uniformly  distributed  sinusoidal  impact  load.  Since 
the  first  two  loadings  are  continuously  applied  over  the  computational  domain  during  the 
analysis  the  converged  transient  solution  is  different  for  each  case.  The  plots  show  the 
differences  between  the  loading  conditions  on  transient  response  effects  for  the  different 
lamination  schemes. 


Figure  13.  Effects  of  plate  theories  on  the  nonlineair  transient  responses  for  the  S3nnmetric 
cross-ply  thick  (a/7i=10)  laminates  with  SSSS  boundary  conditions  and  under  uniformly 

distributed  load  =5x10^). 


Time,  t  (sec) 

(a)  Sym,  cross-ply  laminate  (a/h=10) 


Time,  t(sec) 

(b)  Sym.  cross-ply  laminate(a/h=  100) 


Figure  14(a,b).  Effects  of  loading  conditions  (Impact),  Uniformly  distributed  loading  on  the 
transient  responses  with  SSSS  boundary  conditions  (cross-ply  laminates). 
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Time,  t  (sec)  Time,  t  (sec) 

(c)  Sym.  angle-ply  laminate  (a/h^\0)  (d)  Sym.  angle-ply  laminate  (a/h==100) 

Figure  14(c,d).  Effects  of  loading  conditions  (Impact),  Uniformly  distributed  loading]  on  the 
transient  responses  with  SSSS  boundary  conditions  (angle-ply  laminates). 


CONCLUSIONS 

The  development  of  physically-based  damage  models  and  nonlinear-finite  models  based  on 
the  first-order  shear  deformation  theory  and  the  computer  implementation  of  the  nonlinear 
finite  element  model  were  the  major  accomplishments  of  the  study.  The  integration  df  these 
elements  into  the  layerwise  theory  for  global-local  analysis  and  prediction  of  damage  and 
failures  in  composite  laminates  can  be  carried  out.  The  effects  of  the  lamination  scheme, 
types  of  load,  and  boundary  conditions  on  the  static  and  dynamic  deflection  was 
investigated.  It  is  observed  that  the  effect  of  nonlinearity  on  the  transient  responses  is  to 
decrease  the  amplitude  and  increase  the  frequency  of  nondimensionalized  center 
displacement.  It  so  also  observed  that  the  nonlinear  transient  behaviors  between  TSDT  and 
other  two  theories  (CLPT  and  FSDT)  are  apparent  because  of  the  geometric  nonlinearity. 


RELATIONSHIP  TO  FUTURE  RESEARCH  AND  DEVELOPMENT 

The  integration  of  physically  based  damage  models  that  account  for  micros-scale  as  well  as 
damage  into  robust  computational  models  is  perhaps  the  most  significant  contribution  of 
this  research.  From  the  science  and  technology  points  of  view,  the  results  of  this  research 
will  contribute  to  a  fundamental  understanding  of  the  non-linear  behavior  and  damage 
progression  in  thick  composite  and  sandwich  structures  with  embedded  sensors  and 
actuators.  The  results  will  also  have  a  significant  impact  on  the  design  of  composite 
structures  used  in  lightweight  composite  and  sandwich  structural  components  used  in  army 
and  marine  structures. 
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